Molecular motions in lipid bilayers studied by the neutron backscattering technique 
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We report a high energy-resolution neutron backscattering study to investigate slow motions on 
nanosecond time scales in highly oriented solid supported phospholipid bilayers of the model sys- 
tem DMPC -d54 (deuterated l,2-dimyristoyl-sn-glycero-3-phoshatidylcholine), hydrated with heavy 
water. This technique allows to discriminate the onset of mobility at different length scales for the 
different molecular components, as e.g. the lipid acyl-chains and the hydration water in between the 
membrane stacks, respectively, and provides a benchmark test regarding the feasibility of neutron 
backscattering investigations on these sample systems. We discuss freezing of the lipid acyl-chains, 
as observed by this technique, and observe a second freezing transition which we attribute to the 
hydration water. 



PACS numbers: 87.14.Cc, 87.16.Dg, 83.85.Hf, 83.10.Mj 



I. INTRODUCTION 

Lipid membranes as model systems for more com- 
plex biological membranes^ cannot be understood with- 
out taking into account the structure and dynamics of 
their aqueous environment. The structure and dynamical 
properties of the bound water layers next to the bilayer as 
well as the 'free' (bulk) water further away from the wa- 
ter/lipid interface are of importance in understanding the 
thermal, elastic and transport properties of membranes. 
Furthermore, the interaction between two bilayers is me- 
diated by the hydration watej?2»Milii£. A recent Molecular 
Dynamics (MD) simulation pointed out the importance 
of hydration water dynamics for the understanding of 
the dynamical transition of proteins^. Complementing 
well established structure-function relationships in bio- 
physics, possible dynamics-function relationships remain 
to date much less elucidated. The dynamical properties 
of membrane bound water present an important exam- 
ple in this context. While most spectroscopic techniques 
such as nuclear magnetic resonance (NMR) or dielec- 
tric spectroscopy, are limited to the center of the Bril- 
louin zone and probe macroscopic responses, neutrons 
and within some restrictions also x-rays give unique ac- 
cess to microscopic dynamics on length scales of e.g. 
intermolecular distances. Here, we report on a high 
energy-resolution neutron backscattering study to inves- 
tigate slow motions on nanosecond time scales in highly 
oriented solid supported phospholipid bilayers of the 
model system DMPC -d54 (deuterated 1,2-dimyristoyl- 
sn-glycero-3-phoshatidylcholine), hydrated with heavy 
water. The spectrum of fluctuations in biomimetic and 
biological membranes covers a large range of time and 
length scalesiSiiiiSiiiii^ii&iliiii, ranging from the long 
wavelength undulation and bending modes of the bilayer 
with typical relaxation times of nanoseconds and lateral 
length scales of several hundred lipid molecules to the 
short wavelength density fluctuations in the picosecond 



range on nearest neighbor distances of lipid molecules. 
Local dynamics in lipid bilayers, i.e. dynamics of indi- 
vidual lipid molecules as vibration, rotation, libration 
(hindered rotation) and diffusion, has been investigated 
by incoherent neutron scattering2*i2iiii and nuclear mag- 
netic resonanc o 16 : to determine, e.g. the short wave- 
length translational and rotational diffusion constant. 
Collective undulation modes have been investigated us- 
ing neutron spin-echo spectrometersi^i^ and dynamical 
light scattering^-. Recently, the first coherent inelastic 
scattering experiments in phospholipid bilayers have been 
performed using inelastic x-raji2£ and neutronSi scatter- 
ing techniques to determine the short wavelength dis- 
persion relation. Note that only scattering experiments 
give wave vector resolved access to dynamical properties, 
what is important to associate relaxation times with spe- 
cific motions. 

Information about fluctuations on mesoscopic length 
scales can be inferred from x-ray lineshape analysis 22,23 
in isotropic lipid dispersions. Off-specular x-ray and neu- 
tron reflectivity from aligned phases presents the addi- 
tional advantage that the components of the scattering 
vector Q can be projected onto the symmetry axis of the 
membrane^ 2 ^*^. In both examples, the time-averaged 
elastic scattering is studied, and information on e.g. elas- 
ticity properties and interaction forces can be obtained. 

Contrarily, dynamical properties like transport coeffi- 
cients can be inferred from direct measurements of dy- 
namical modes. Again, these measurements are prefer- 
ably carried out in aligned phases to preserve the unique 
identification of modes on the basis of the parallel and 
perpendicular components Q r and Q z of the scattering 
vector Q, as pictured in Fig. 2] Recentlj^i, we have 
demonstrated that collective supra-molecular dynamics 
of planar lipid bilayers, notably the dispersion relation 
of density modes in the lipid acyl-chains, can be stud- 
ied using the so-called three-axis neutron spectroscopy 
technique giving access to an energy resolution of up to 
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Monochromator 

FIG. 1: (Color online) Schematic of the scattering geometry 
(not drawn to scale) . The sample was mounted with the mem- 
brane planes vertical, i.e. perpendicular to the horizontal scat- 
tering plane of the backscattering spectrometer. The mem- 
brane planes were oriented at well-defined angles with respect 
to the incident beam. The inter acyl-chain correlation peak 
in the plane of the membranes is located at Q r = lA2 A -1 . 
Spatially arranged analyzers at a sample-to-crystal distance 
of approx. 1.5 m allow to separately but simultaneously probe 
the dynamics on different length scales. The bulk heavy wa- 
ter correlation peak, which occurs at Q r =2 A -1 , is indicated 
by the dotted line. 



about 300 /zeV . The low scattering volume of quasi two- 
dimensional planar membranes and the small inelastic 
signal pose particular experimental challenges. Sample 
preparation and experimental set-ups have to be adapted 
for inelastic experiments. In Ref. [2l| the scattering vol- 
ume restriction was overcome by stacking several thou- 
sand highly aligned bilayers. These experiments shall 
be complemented in the present work by /zeV energy 
resolved spectra, achieved by the neutron backscatter- 
ing technique. We thereby gain access to the low ener- 
getic density fluctuations corresponding to slow motions 
on nanosecond time scales. By analyzing the respective 
Q dependence, we simultaneously probe internal length 
scales from 3 to 18 A and investigate freezing of the lipid 
acyl-chains and of the hydration water, i.e. the water 
layer in between the stacked membranes. The paper is 
organized as follows: In the next section, we briefly dis- 
cuss hydration water in stacked planar membranes based 
on literature. Sect. Hill gives details of the experimental 
set-ups used for neutron backscattering and diffraction 
on three-axis spectrometers. The corresponding results 
are presented in Sects. IIII Al and IIII Bl and discussed in 
Sect. El 



II. HYDRATION WATER IN MULTILAMELLAR 
ORIENTED MEMBRANE STACKS 

Freezing of water in lamellar structures (DEPE, 
DEPC, DOPC) was first reported by Gleeson et al. 26 
from x-ray diffraction. Konig et al& have investigated 
the molecular dynamics of water in DPPC multilayers by 
quasielastic neutron scattering (QENS) and nuclear mag- 
netic resonance (NMR). Neutron diffraction and QENS 
experiments on the interaction of hydration water and 
purple membranes have been reported by Lechner and 
Fitter and coworkers2L2&. Additionally, Marrink et al. 
have investigated the ordering of membrane water in 
DPPC by Molecular Dynamics simulations (MD) 29 . 

The combination of experimental results and conclu- 
sions of the complementary techniques and experiments 
and the outcome of the simulations gives a plausible pic- 
ture of the behavior of hydration water in multilamel- 
lar oriented membrane stacks. The basic scenario is the 
following 2 -: The lamellar spacing d z proves to be a very 
fine measure for the water content in between the mem- 
brane stacks as it is sensitive to the number of water lay- 
ers (hydration shells) between two bilayers (which have 
a thickness of about 3.5-4 A). The membrane hydration 
water can be cooled well below the normal freezing point 
of (bulk) water. d z shows a hysteresis in heating and cool- 
ing cycles. When on cooling the temperature falls below 
the freezing temperature of bulk water and the hydration 
water becomes supercooled, d z starts to drop indicating 
a decreasing number of water layers between the stacks. 
Ice Bragg peaks appear with an intensity increase propor- 
tional to the decrease of d z . The authors conclude that 
the hydration water migrates out of the membrane stacks 
when freezing and condenses as ordinary ice outside the 
lamellar structure and the remaining hydration water is 
in equilibrium with bulk ice. There is no crystalline ice 
between lipid bilayers in an ordered multilamellar phase. 

From QENS measurements at different degrees of hy- 
dration and orientations (with the scattering vector Q 
parallel and perpendicular to the plane of the bilay- 
ers) , the different water dynamics and anisotropy can be 
extracted 9 . At low hydration of the bilayers, rotation of 
water molecules is the dominant motion without signs of 
translational diffusion. At high hydration a translational 
diffusion has to be allowed to fit the QENS spectra. The 
rotational degree of freedom corresponds to that of free 
bulk water, while the translational diffusion appears to 
be hindered, however, as compared to bulk water with 
a slightly but significantly smaller diffusion constant D. 
The water dynamics is isotropic within the length scales 
probed of about 10 A. The first hydration layers around 
the lipid head-groups are tightly bound and show rota- 
tional dynamics, only. Additional layers participate in 
the translational dynamics and behave like quasifree bulk 
water. The authors estimate that about 40% of the water 
in a hydrated sample is tightly bound. 

The MD simulations show that the water structure is 
highly perturbed by the presence of the two membrane 
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surfaces^. The analysis of orientational polarization, hy- 
drogen bonding and diffusion rates of the water molecules 
between the membranes reveals a strong perturbing ef- 
fect, which decays smoothly towards the middle of the 
water layer. Bulk-like water is only found at a distance 
of 10 A away from the interfacial plane. But even there, 
the hydrogen bonding pattern as well as diffusion rates 
show small but significant deviations, in agreement with 
the experiments. 



III. EXPERIMENTAL 

Deuterated DMPC -d54 (deuterated 1,2-dimyristoyl- 
sn-glycero-3-phoshatidylcholine) was obtained from 
Avanti Polar Lipids. Highly oriented membrane stacks 
were prepared by spreading a solution of typically 
25 mg/ml lipid in trifluoroethylene/chloroform (1:1) on 
2" silicon wafers, followed by subsequent drying in vac- 
uum and hydration from D2O vapor—. Fifteen such 
wafers separated by small air gaps were combined and 
aligned with respect to each other to create a 'sandwich 
sample' consisting of several thousands of highly oriented 
lipid bilayers (total mosaicity of about 0.6°), with a total 
mass of about 400 mg of deuterated DMPC. 

The experiment was carried out at the cold neu- 
tron backscattering spectrometer ENIGMA at the Institut 
Laue-Langevin (ILL) in its standard setup with Si(lll) 
monochromator and analyzer crystals corresponding to 
an incident and analyzed neutron energy of 2.08 meV 
(A=6.27 A) (as sketched in Fig^). Neutron backscat- 
tering was used to obtain wave vector-resolved dynam- 
ical information with a high resolution in energy trans- 
fer of about 0.9 /ieV FWHM. The high energy resolu- 
tion in neutron backscattering results from the use of 
90° Bragg angles at both the monochromator and ana- 
lyzer crystals^. Thus in principle the precision in the 
determination of the neutron wavelength is only limited 
by the Darwin width of the crystals used™. In the exist- 
ing experimental implementations of the backscattering 
geometry, contributions from the beam divergence and 
deviations from exact 90° Bragg angles slighly deterio- 
rate this precision. As is indicated in Fig^ in the IN10 
spectrometer a small deviation from exact backscatter- 
ing occurs at the monochromator (less than 1°) to al- 
low for the returning neutrons to hit the deflector and 
subsequently illuminate the sample inside the secondary 
spectrometer. The scattered neutrons returning from the 
analyzers are detected in exact backscattering and thus 
may travel twice through the sample. This does not pose 
a problem since the total scattering probability by the 
sample is always kept below 10%. The incident beam is 
adequately pulsed, and therefore directly scattered neu- 
trons can be discriminated from analyzed neutrons by 
their flight time. 

The incident beam at the sample in our experiment 
was about 3 cm wide and 4.5 cm high with a diver- 
gence corresponding to the critical angle of Ni 58 . The 



sample was mounted in a hermetically sealed aluminium 
container within a cryostat and hydrated from D2O va- 
por. Saturation of the vapor in the voids around the 
lipids was assured by placing a piece of pulp soaked in 
D2O within the sealed sample container. The pulp was 
shielded by Cadmium to exclude any parasitic contribu- 
tion to the scattering. The hydration was not controlled 
but we allowed the sample to equilibrate for 10 h at room 
temperature before the measurements. We note that the 
absolute level of hydration of the sample is therefore not 
precisely knowni^, but the temperature of the main tran- 
sition agrees quite well with literature values. We argue 
that the question of hydration might become more im- 
portant when discussing relaxation times and diffusion 
constants while in this study we mainly concentrate on 
phase transitions and elastic scattering. 

Using the IN10 spectrometer, two types of measure- 
ments could be performed: Firstly, fixed energy-window 
scans centered at zero energy transfer (FEW-scans), have 
been recorded as a function of the sample temperature. 
This scattering intensity can be regarded as purely elas- 
tic, within the excellent instrumental resolution render- 
ing motions faster than approx. 4 ns already visible. 
Changes in this intensity may arise from either structural 
or dynamical changes, i.e. shifts of correlation peaks or 
freezing of dynamical modes. In connection with diffrac- 
tion data, structural changes may be separated from dy- 
namical effects. Thus, from FEW-scans, information on 
the onset and type of molecular mobility in the sample 
can be inferred and glass or melting transitions clearly 
identified and assigned to corresponding length scales by 
analyzing the corresponding Q-dependence. The second 
type of measurement was performed by Doppler-shifting 
the incident neutron energy through an adequate move- 
ment of the monochromator crystal. These energy scans 
correspond to a time range of motion in the sample on 
the order of 10~ 9 s < t < 10~ 8 s. 

In view of its entangled geometry, the application of 
the neutron backscattering technique for probing dynam- 
ics at interfaces is challenging. Nevertheless, we demon- 
strate here the feasibility of backscattering on the lipid 
membrane sandwich samples. Analogously to the three- 
axis experiments, the samples have been oriented in the 
spectrometer to measure at wave vector transfers parallel 
and perpendicular to the lipid membrane plane, respec- 
tively. The IN 10 analyzers cover an angular range of 
approximately 20° each, resulting in a rather poor ir- 
resolution (which does not allow to measure diffraction 
data to determine, e.g. the lamellar spacing d z ), but en- 
hanced sensitivity for even very small inelastic signals. 
Note that we are therefore not sensitive to shifts of cor- 
relation peaks within this coarse Q resolution and cor- 
responding structural changes of the bilayers. The six 
discrete detector tubes (D1-D6) of IN10 cover a total Q- 
range of 0.34 A^ 1 up to 1.91 A~\ namely 
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-l 


1.59 A- 


1 <D5<1.78 A^ 1 


D6 


1.85 A" 


-l 


1.78 A~ 


1 <D6<1.91 A -1 . 



The onset of molecular mobility can thus be mea- 
sured separately but simultaneously for different molec- 
ular components at different Q values and length scales. 
The lipid acyl-chain correlation peak that occurs at 
Q r ~1.42 A -1 was mainly detected in D4, as pictured 
in Fig.n Note that the diffraction angle of (heavy) bulk 
water (Q=2 A -1 ) is not directly accessible on IN10 for 
geometrical reasons. The detector at highest angle (D6), 
which is centered at Q=1.85 A -1 , covers a Q range of ap- 
prox. 1.78 A -1 < Q <1.91 A -1 and measures in the tails 
of the broad (heavy) water correlation peak. D5 and 
D6 also cover the range of the (hexagonal, PQ^/mmc) 
ice Bragg peaks (100), (002) and (101) at Q values of 
1.605 A" 1 , 1.706 A" 1 and 1.817 A" 1 , respectively. We 
have measured all scans (elastic and inelastic) also with 
the membranes rotated by 90°, hence with Q perpendic- 
ular to the membrane surface (Q z ), as a reference. 

Complementary diffraction data have been measured 
on the cold and thermal three-axis spectrometers IN12 
and IN3. IN3 was equipped with its multi analyzer 
detector— to map out large areas of reciprocal space si- 
multaneously, with all analyzers aligned to a fixed final 
energy of the scattered neutrons of i?/=31 meV. IN12 
was used with fixed kf =1.75 A -1 resulting in a Q reso- 
lution of AQ=0.02 A -1 . For the three-axis experiments, 
the sample was held in a dedicated humidity chamber 
and - as for the backscattering measurements - hydrated 
from the vapor phase. The use of a spectrometer to mea- 
sure diffraction (energy transfer Auj = 0) has a partic- 
ular advantage when working with partially deuterated 
and/or biological samples containing C, O, N, H (or D). 
As the downstream analyzer cuts out the (within the in- 
strumental resolution) elastically scattered neutrons, the 
quasielastic contribution to the background arising from 
incoherent scattering is reduced and the signal to noise 
ratio drastically improved. 



A. Elastic neutron diffraction 

Figure [21 shows Q r scans of DMPC, measured on IN12 
and IN3 at T=18°C, for a dry (ambient humidity) and 
a hydrated (d 2 =58.5 A -1 ) DMPC sample, respectively. 
The measured Q r rang e of 0.65 A" 1 < Q r < 2.25 A" 1 
covers the correlation peaks of the lipid acyl-chains (Q r — 
1.42 A -1 ) and heavy bulk water (Q r ~ 2 A -1 ), which 
originate from the average nearest neighbor distances of 
acyl-chains and water molecules, respectively. The addi- 
tional peak observed at around Q r ~ 0.7 A -1 could be 
a contribution of the packing of the lipid head groups in 
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FIG. 2: (Color online) Diffraction data of a dry (ambient 
humidity) and a hydrated sample (T=18°C) taken on IN12 
and IN3, respectively. The scans cover the Q range of the 
correlation peaks of the lipid acyl-chains and of water. The Q 
coverage of the six discrete IN10 detectors D1-D6 is indicated 
in the Figure. (Dotted lines mark positions of Aluminum 
powder lines.) 

these gel or sub-gel phases, which is usually not seen with 
x-rays. But the exact assignment will need a dedicated 
study. As is well known, the Q r position of the chain 
correlation maximum shifts to smaller values upon hydra- 
tion, reflecting the corresponding increase of the nearest 
neighbor distance. An additional broad peak appears at 
around Q r ~2 A" 1 (FWHM 0.4 A" 1 ) and is attributed 
to the water layer between the membranes. Figure |3| 
shows a mapping of the elastic scattering intensity in the 
Qr/Qz plane in the hydrated gel phase, at T=18°C. The 
intensity distribution reflects the quasi two-dimensional 
liquid structure, as discussed in Ref. l35l The elongated 
rod-like peak can be understood by analogy with the well 
known crystal truncation rods in surface crystallography. 
This smearing of the intensity in the perpendicular direc- 
tion is also well known from scattering of monolayers at 
the air water interface. Contrarily, the weak intensity at 
the water position (Q r — 2 A -1 ) is truly isotropic and 
distributed over a Debye-Scherrer ring. Consequently, we 
expect a strong anisotropy in the backscattering exper- 
iment at the Q positions of the lipid acyl-chains, while 
scattering at the water position should be isotropic. 

B. Backscattering 

The intensity measured in the FEW scans is essentially 
equivalent to the 'elastic structure factor' S{Q,U) = 0), 
except for minor temperature independent corrections, 
e.g. due to self-absorption or geometrical effects. In con- 
trast to the x-ray or neutron elastic scattering factor, 
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FIG. 3: (Color online) Q r /Q z mapping of the diffraction sig- 
nal from the hydrated sample (gel phase, T=18°C). The cor- 
responding Q z = section is shown in Fig.[2](IN3 multidetec- 
tor data). The correlation peak of the acyl-chains exhibits the 
shape of a slightly bent Bragg rod, reflecting the quasi two- 
dimensional liquid-like short range correlations of the acyl 
chain positions. 



this is the 'true' elastic scattering component which is 
measured with the highest possible energy resolution and 
minimal contribution from quasi-elastic or inelastic scat- 
tering. Changes in this intensity are basically sensitive 
to dynamical changes, since the excitation of dynamical 
modes will shift scattering intensity from zero to finite en- 
ergy transfer, out of the tight fixed energy window. Note 
that FEW data may include contributions from coherent 
and incoherent scattering, i.e. from pair- and autocorre- 
lated scattering (collective and local modes). 

Figure 0] shows FEW scans for all detectors in a tem- 
perature range of 200-315 K to map out the transition 
of the different molecular components from immobile to 
mobile as a function of temperature for the two set-ups, 
with (a) the scattering vector Q placed in the plane of the 
membranes and (b) perpendicular to the bilayers. Both 
measurements have been taken while cooling with a con- 
stant rate of 0.1 K/min and counting times of 5 min per 
point. Note that we are therefore not sensitive to possible 
temperature shifts smaller than the temperature resolu- 
tion of about 0.5 K. The data are shown as raw data 
without any correction, solely normalized to the incident 
beam monitor. The contributions by the Aluminum can 
and cryostat and Silicon wafers are in any case tempera- 
ture independent, as is the individual detector efficiency 
and effects induced by the sample geometry. Note that 
the wavelength used is too long to excite the Aluminum 
Bragg peaks. The in-plane component (Q r ) in Fig. 0] 
(a) shows a distinct step in D4 at Tf =293 K, which 
is the temperature of the main transition from the P/3' 
(gel) into the fluid La phase of the deuterated DMPC 
bilayers. The transition occurs some degrees lower in the 



deuterated compound (at 20.15°C=293.3 K 36 ) as com- 
pared to protonated DMPC and is detected in D4, where 
the maximum of the static structure factor of the lipid- 
chains (the acyl-chains correlation peak) occurs. A sec- 
ond, kink-like anomaly appears in detectors D5 and D6 
at about Tf w =271 K. In the perpendicular direction (Q z 
shown in Fig.0](b)), D5 and D6 display a similar temper- 
ature behavior as for Q r . Thus, the scattering at the cor- 
responding Q position can be concluded to be isotropic. 
Contrarily, the two D4 curves along Q z and along Q r , 
respectively, exhibit different functional forms. In par- 
ticular the step-like anomaly at Tf =293 K is lacking in 
the Q z curve, indicating an anisotropic scattering at this 
Q. This strong anisotropy of the acyl-chain correlations 
has been suggested by the elastic intensity distribution in 
Fig. |21 It is therefore not astonishing that temperature 
induced changes in the dynamics (within the experimen- 
tal window of time and length scales) measured at the Q 
position of the acyl chain peak is also anisotropic. 

To quantify the anisotropy of the elastic scattering, we 
consider the intensity along Q z as 'background' (Fig. 0] 
(b)) and subtract these scans from the scans in Fig. 0] 
(a) . Thus in Fig. [3] we plot the monitor normalized in- 
tensity differences I n orm(Qr) — InormiQz)- In this repre- 
sentation, isotropic contributions to the signal cancel out, 
and anisotropic contributions become more pronounced 
due to the background subtraction. Note that this in- 
tensity difference becomes negative if the elastic contri- 
bution along Q z is stronger than along Q r and can no 
longer to be considered as a background. This is the 
case for the small angle detectors D1-D3, which are ex- 
pected to pick up signal from the lamellar Bragg sheet 
tails (non-specular reflectivity). As shown in Fig. [5] two 
well separated freezing steps ('immobile' within the res- 
olution window), of the lipid acyl-chains at Tf =293 K 
and at about 271 K are clearly visible in detectors D4, 
D5 and D6. Note that D5 obviously still shows traces of 
the freezing step in the lipid acyl-chains. 

Figure shows the important temperature range 
around the main phase transition for the D4 detector, 
which is dominated by the lipid acyl-chain correlation 
signal. There is a sharp freezing step at T/=293 K, 
corresponding to the temperature of the gel-fluid phase 
transition in deuterated DMPC. Above Tf up to about 
302 K, we observe a regime with a smaller slope which 
coincides with the temperature range of pseudocritical 
(anomalous) swelling in DMPC, as it will be discussed 
below. 

Fig.[7|shows representative energy transfer scans of the 
deuterated DMPC membrane stack. The data have been 
taken at three different temperatures, at T=250 K, 290 K 
and 300 K with a typical counting time of about 9 hours 
per temperature. For these scans, the detected neutrons 
are sorted into discrete channels which are evenly spaced 
in the measured voltage signal arising from an induction 
coil at the mechanical Doppler drive. The energy transfer 
is then calculated from this voltage signal by using the 
known frequency of the sinosoidal Doppler monochroma- 
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FIG. 4: (Color online) Fixed energy-window scans with the scattering vector Q placed in the plane of the membranes (a) and 
perpendicular to the bilayers (b). The curves are shifted along the intensity axis for clarity. The drawings show the orientation 
of the scattering vector Q with respect to the membrane stacks. (Counting is normalized to monitor) 



tor movement and the constant channel width in terms 
of voltage. Therefore, we note that a systematic error 
may not be completely ruled out for the energy scale of 
the data. From the data, the onset of molecular mobility 
on the corresponding length scales can be inferred. The 
presence of an elastic peak in the spectra points to static 
order over the experimental time window. Note that a 
fluid system does not show elastically scattered inten- 
sity, i.e. no order at infinitely long time scales. Even 
within the very limited statistics, the different dynam- 
ics are clearly visible: While the lipid acyl-chains melt 
between 290 and 300 K, melting at the water position 
already occurs between 250 and 290 K. 



IV. DISCUSSION 

The backscattering detectors cover the Q range of the 
correlation peaks of the lipid acyl-chains, bulk water and 
ice. The diffraction data in Figs. |21 and |3] show that 
the corresponding correlations lead to signals in the neu- 
tron scattering experiments and allow to study simul- 
taneously lipid and water dynamics in the backscatter- 
ing experiment at the positions of the corresponding cor- 
relation peaks. In particular, this study allows to dis- 
criminate the immobile-to-mobile transitions of the lipid 
acyl-chains and the water. The dynamics of the mem- 
brane associated water is relevant for biological mem- 
brane function in the physiological temperature range. 
The neutron backscattering technique can contribute to 
a quantitative understanding of the water dynamics in 
membranes. Structural and dynamical changes can lead 



to an increase or loss of elastically scattered intensity in 
the backscattering experiment. A structural transition or 
change might lead to a shift of correlation peaks, which 
then move out of the respective detector. Slowing down 
(freezing with respect to the instrumental resolution and 
the accessible time window) of dynamical degrees of free- 
dom gives a quasi-static contribution and leads to an in- 
crease in the elastically scattered intensity. Because of 
the relaxed Q resolution, the backscattering experiment 
does not allow to probe small Q-shifts of the acyl chain 
correlation peak^i and we therefore argue that we are 
basically sensitive to dynamical changes. 

There is a sharp freezing step in the FEW scans in 
the detector located at Q=1.42 A" 1 , indicating a first 
order transition at 7/ =293 K. The onset of mobility in 
the acyl-chains of the lipid membranes probably devi- 
ates from a simple melting transition, as suggested by 
the slow decrease of the elastically scattered intensity 
with temperature in the range 250 K<T<290 K (see 
Figs. 1415( 1 . The melting (or freezing) of the lipid acyl- 
chains in DMPC and other lipids is indeed of first or- 
der, but is known to show a pseudocritical swelling, i.e. 
a continuous change of the lamellar d^-spacing above Tf 
(see Refs. I37ll38ll39ll40r) and also of the inter acyl-chain 
correlation peakSi. The backscattering data probes the 
temperature induced changes of dynamical modes with 
ultra high energy but very modest Q resolution. Here 
we observe a region above Tf of Tf < T < Tf + 9 K, 
with a smaller slope in the anisotropic elastic scattering 
(Fig. [fJJl which coincides well with the region of pseud- 
ocritical swelling. There is an ongoing discussion about 
(pseudo)criticality in lipid bilayers and the different con- 
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FIG. 5: (Color online) Anisotropy of the elastic scattering sig- 
nal (see text for explanation). The mobile-immobile transi- 
tion (within the resolution window) is clearly different for the 
'lipid chains detector' and the 'water peak detectors'. Charac- 
teristic steps indicating the excitation of dynamical modes (or 
freezing of modes upon cooling) are observed at Tf =293 K (at- 
tributed to lipid acyl chains) and at J 1 /™ =271 K (attributed 
to the membrane water). Solid lines are guides to the eye. 
(Counting is normalized to monitor). 



tributions of water layer and lipids to the anomalous ex- 
pansion. Pabst et a0~- report from x-ray diffraction on 
DMPC bilayers that the anomalous swelling is essentially 
the result of an expansion of the water layer and caused 
by increased fluctuations. They find a softening of the 
bilayer in the vicinity of Tf, i.e. a (small) increase of the 
Caille fluctuation parameter rji. The present backscat- 
tering study allows to assign the fluctuations to specific 
length scales. The critical swelling regime is thus accom- 
panied by significant changes in the local and collective 
in-plane dynamics (diffusion, density fluctuations, undu- 
lations) of the lipid acyl-chains in the whole range of 
pseudocriticality. 

A quantitative analysis of the step height in Fig. [B] 
shows that about 85% of the 2D dynamical modes in the 
plane of the membranes freeze at the main transition at 
Tf. 15% freeze already in the region of critical swelling 
of the bilayers which may be accompanied by coexisting 
gel and fluid microdomains, as discussed in Ref. |^| and 
also reported from NMR experiments 9 . 

We attribute the second freezing transition at 
Tf w — 271 K, which is observed in D5 and D6, to freezing 
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FIG. 6: (Color online) Freezing of dynamical modes in the 
lipid acyl-chains, observed in detector D4 (Q = 1.42 A -1 ). 
There is a sharp intensity step at Tf =293 K, corresponding to 
the temperature of the gel-fluid phase transition in deuterated 
DMPC. Above Tf up to about 302 K there is a temperature 
region characterized by different, less steep slope (marked by 
the shaded area). This temperature regime is also denoted as 
the anomalous swelling regime, since the lamellar periodicity 
d z is known to show a continuous swelling in this range. 



of the membrane hydration water. Upon cooling, mo- 
bile bulk like water in the center of the water layer may 
leave the layer in between the membranes and condense 
as polycrystalline ice outside of the lamellar structure. 
Freezing of this water leads to disappearing of the respec- 
tive dynamic modes and is responsible for the sharp freez- 
ing step in Fig. [3] The polycrystalline bulk ice then gives 
rise to an elastic contribution in detectors D5 and D6, 
as it is observed in Fig. Unlike the lipid acyl-chains, 
the water shows a more gradual freezing. Freezing of the 
hydration water is lowered by about six degrees, as com- 
pared to pure heavy water ( 2 H 2 0) at 276.97 Ki 2 . The 
reduced number of hydrogen bonds for the interstitial 
water, which was found in MD simulations 2 ^, might be 
responsible for the supercooling. An interesting result of 
the present experiment is the fact that after subtraction 
of the FEW scans along the two principal axis of lamellar 
symmetry, there is a residual excess scattering contribu- 
tion also in the so-called 'water peak' detectors. In other 
words, a fraction of the water exhibits anisotropic scat- 
tering, and the freezing step of this component can be 
probed. 

Below Tf W) all detectors in Fig. 0] show an increase 
in the elastically scattered intensity, ft seems that the 
membranes become 'stiffer' on all observed length scales 
upon cooling. This reduction of dynamical modes may 
also be linked to the decrease of the inter-bilayer distance 
when water freezes out. 
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FIG. 7: (Color online) Energy scans at temperatures 
T=250 K, 290 K and 300 K for the Q-values 1.42 A" 1 (lipid 
acyl-chain correlation peak) and Q=1.85 A -1 (hydration wa- 
ter correlation peak). At 290 K, the water is already 'mobile' 
within the experimental energy resolution whereas the lipid 
acyl-chains are still frozen. (Counting is normalized to moni- 
tor) 



V. CONCLUSION AND OUTLOOK 

We present the first neutron backscattering study to 
investigate molecular mobility of the different molecular 
components in a phospholipid membrane stack. By com- 
bining backscattering and diffraction data, we are able to 
separate pure structural from dynamical changes. The 
main transition from the gel to the fluid La phase is con- 
nected with freezing of collective or diffusive motions of 
the acyl-chains. In the temperature range of the so-called 
pseudocritical (anomalous) swelling we find a continuous 
tail in the FEW scans at the high temperature side of the 
freezing step, presenting a deviation from the simplest 
scenario of a first order freezing transition. A second 



transition is observed at Tf w ~271 K which we attribute 
to water molecules in the layer between the membrane 
stacks. 

The experiment also provides a benchmark test for 
inelastic neutron studies of planar, i.e. quasi two- 
dimensional membrane model systems. Neutron beams 
for high resolution experiments are weak even at the most 
intense neutron sources, and this restriction of the low 
scattering intensity has been successfully circumvented 
by stacking the membranes as described above. The 
data are, however, still incomplete. In particular, the 
low count rate in the individual detector tubes does not 
allow for sufficient statistics (Fig. UJ. Therefore, future 
experiments will be aimed at testing whether a quasi- 
elastic broadening indicating slow motion on nanosecond 
time scales occurs in the planar lipid membranes. In 
combination with the three-axis technique, where col- 
lective short wavelength fluctuations in the picosecond 
range are probed, the backscattering technique could 
thus give access to relaxations on the same length, but 
in the nanosecond time range. We also note that in the 
present study both the lipid acyl-chains and the water 
were deuterated. The experiment is therefore predom- 
inantly sensitive to detecting collective motions arising 
from coherent scattering. Further experiments will in- 
clude a selective deuteration of the acyl-chains and the 
membrane water, respectively, to mask different types 
of mobility. In addition, structural information will be 
obtained by using diffraction detectors in parallel with 
the energy-discriminating backscattering detectors. The 
position of the reflectivity Bragg peaks and the lipid acyl- 
chain correlation peak will allow simultaneous detection 
of the lamellar spacing d z and thickness of the water layer 
d w , as well as the ordering of lipids in the plane of the 
membranes. 
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